Correlated metallic two particle bound states in quasiperiodic chains 
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Single particle states in a chain with quasiperiodic potential show a metal-insulator transition 
upon the change of the potential strength. We consider two particles with local interaction in the 
single particle insulating regime. The two particle states change from being localized to delocalized 
upon an increase of the interaction strength to a nonperturbative finite value. At even larger 
interaction strength the states become localized again. This transition of two particle bound states 
into a correlated metal is due to a resonant mixing of the noninteracting two particle eigenstates. 
In the discovered correlated metal states two particles move coherently together through the whole 
chain, therefore contributing to a finite conductivity. 
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The interplay of localization and many body interac- 
tions has been a highly active research topic since the 
discovery of Anderson localization pQ. While the direct 
theoretical study of such systems is quite complex and re- 
cently hotly debated (see e.g. [2 ), another route is taken 
by studying few interacting particles. The potential ap- 
plicability to recent experimental activities with ultracold 
atoms in optical lattices [3] increased the interest in cor- 
responding theoretical studies. Localization of single par- 
ticle states in a one-dimensional lattice can be achieved 
using a lattice potential which is biased by an external 
dc field leading to Wannier- Stark localization and Bloch 
oscillations [5], random leading to Anderson localization 
PP, and quasiperiodic leading to Aubry- Andre localiza- 
tion [4]. Two interacting particles (TIP) in a dc field lead 
to no substantial change of the localization length due to 
the Stark ladder structure of the single particle eigenener- 
gies (disregarding some resonant tunneling events which 
however are also suppressed for large distances) [6]. TIP 
in the random case do lead to an increase of the localiza- 
tion length with some controversial discussions about the 
quantitative outcome [7-9 , but do not yield complete de- 
localization, For TIP in the quasiperiodic potential case 
few numerical results give varying predictions from incre- 
mental increase of localization length to opposite reports 
of decrease of localization length in the insulating regime 



In this work we consider the TIP problem in a 
quasiperiodic chain at finite (nonperturbative) strength 
of interaction, deep in the single particle insulating 
regime. We observe a complete derealization of cer- 
tain two particle bound states, in which both particles 
are keeping a relative distance less than the one par- 
ticle localization length. The interaction renormalizes 
eigenenergies of different classes of localized TIP eigen- 
states with different strength, and leads to a resonant 
overlap of these energies in a certain range of the in- 
teraction constant. In this nonperturbative window the 
overlap between these groups of eigenstates leads finally 



to a complete derealization. This results in the novel 
observation of a correlated metal state build from only 
two interacting particles. 

We study the TIP in the framework of the Hubbard 
model with Hamiltonian 
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where b^ and bj are creation and annihilation operators 
of indistinguishable bosons at lattice site j, and U mea- 
sures the onsite interaction strength between the bosons. 
The potential 
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controls the single particle problem via its strength A. 
The incommensurability parameter a = {y/h — l)/2 
(golden mean). j3 is an arbitrary phase which controls 
wave packet dynamics and is irrelevant for the properties 
of extended metallic eigenstates. 

For a single particle the interaction term does not con- 
tribute. Using the basis \j) = |0) the eigenstates \v) 
of Hamiltonian (HI) with eigenvalues \ v are computed us- 
ing \v) = E^&flO) , U\v) = K\u). A) = 0» are 
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the eigenvectors. As was first shown by Aubry and An- 
dre [4] (see also [13j E]), for A < 2 all eigenstates are 
delocalized (metallic phase) and for A > 2 all states are 
localized (insulating phase), with a localization length 
£i = l/ln(A/2). Consequently the probability distribu- 
tion function (PDF) = (q\b+bi\q) ~ exp [-2|Z|/fi] 
is exponentially localized in the insulating regime. The 
bounded eigenvalue spectrum shows fractal properties 
(see e.g. [15]. At the value A = 2.5 (which is the main ref- 
erence parameter for computational studies in this work) 
the spectrum has three main minibands called SP1, SP2, 
and SP3 with centers around A = — 2.5, 0, 2.5 and ordered 
with increasing energy. 

For two particles we expand the eigenstates \q) of the 
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TIP problem in the local basis 
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(l,m\q) are the normalized eigenvectors 
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with Z < m. Thus [£j^J 2 is the probability to find two 



particles on the sites I and m. We also compute the PDF 
of the eigenmodes \q) as 
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and their participation numbers P q = l/J^f" • 
We note that at the noninteracting limit U = and 
A = 2.5 the spectrum decomposes into five minibands 
TP1,TP2,TP3,TP4,TP5 (see Fig. [TJbelow). This follows 
from the three miniband structure of the single particle 
spectrum, such that e.g. TP1 is formed from two single 
particle states SP1: TP1 = (SP1 x SP1). Analogously, 
TP2 = (SP1 x SP2), TP3 = (SP2 x SP2) U (SP1 x 
SP3), TP 4, = (SP2xSP3), TP5 = (SP3xSP3), where 
TPp = (SPmxSPn) corresponds to a two-particle prod- 
uct state TPp with one of the particles being in a single- 
particle eigenstate in miniband SPm and another one in 
SPn. 

To check for the selective character of a derealiza- 
tion effect we present participation numbers of all eigen- 
states versus their eigenenergies for the particular choice 
of interaction constant U = 7.9 (Fig. []Ja)). The five 
miniband structure of the noninteracting case is clearly 
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FIG. 1: (a) Participation number of interacting particle eigen- 
states versus their eigenenergies, (b) and (c) present typical 
localized and delocalized eigenfunctions \C[ q l\ which corre- 
spond to the background and peaks in graph (a), respectively. 
In (d) we plot PDFs of localized (red) and delocalized (blue) 
eigenmodes. Here U = 7.9 and A = 2.5. 
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FIG. 2: Dependence of the largest participation number (of 
the mostly delocalized eigenstate) on U for different lattice 
sizes N = 100, 150, 200, 250 (from bottom to top). Here A = 
2.5. 



seen. This follows from the fact, that most states cor- 
respond to two particles separated by a distance larger 
than the localization length, and therefore these states do 
not change when the interaction is increased. It is these 
states which also stay at a small value of the participa- 
tion number of the order of P ~ 5. However in the mini- 
bands TP2 and TP4 we observe candidates for metallic 
delocalized states with 50 < P < 100 being one order 
of magnitude larger. Two characteristic eigenvectors for 
localized and delocalized states are shown in Fig|I'b,c), 
and their corresponding PDFs are plotted in Fig lid). 
The metallic state in Fig|TJb,d) is indeed occupying the 
whole system. For this state both particles stay close 
to each other, forming a diagonal structure in Figjljb). 
Therefore we coin these correlated states metallic two 
particle bound states. 

Next we examine the characteristics of eigenstates for 
A = 2.5 and various systems sizes and different values 
of U. We start with N = 100. For each value of U we 
find the state with the largest participation number and 
plot this number versus U in Figj2] While for U < 4 and 
U > 10 the number stays around 10-20, three humps up 
to values of 40 are noted for 4 < U < 10. Since the 
spatial size of an eigenvector is roughly 2-3 times larger 
than its participation number, we conclude that in the 
hump regions the analyzed states extend over the whole 
system. We increase the system size to N = 150, 200, 250 
and repeat the above analysis. We observe that the hump 
heights grow linearly with the system size, indicating that 
for all system sizes the most delocalized eigenstates in 
the hump regions are spreading over the whole system. 
This is a clear evidence of complete derealization of some 
eigenvectors in the mentioned parameter region. 

Using exact diagonalization we are restricted to lattice 
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FIG. 3: (Color online) Time and space dependence of the 
PDF of an initial state with two particles at adjacent sites for 
A = 2.5 (color maps log 10 PDF). Left panel: isolator phase, 
U = 2, right panel: correlated metal phase, U = 4.5. 
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sizes N < 300. In order to push the limits, we compute 
the evolution of the Schrodinger equation in real time 
without diagonalization, starting with two particles lo- 
cated on adjacent sites. When and if the extended states 
exist, such initial conditions must excite them, making 
a part of a two-particle wave packet propagate. Typical 
PDFs as a function of space and time are shown in FigJ3] 
We find ballistic spreading over the whole system with 
N = 2500 sites for U = 4.5, and complete localization 
for U = 2. This extends the evidence for complete der- 
ealization of the initial state into extended eigenstates of 
two interacting particles. 

We perform a scan in the parameter space {[/, A} in 
order to identify the region of correlated metallic two 
particle bound states. We choose a system size N = 610, 
and place two particles at adjacent sites in the center of 
the chain. The Schrodinger equation is evolved up to 
time t = 1.5 x 10 4 . The square root second moment of 
the wave packet PDF is then measured for 60 different 
original particle positions and the outcome for the fastest 
growing realization is plotted in color code as a function 
of both A and U in FigjiJ For A < 2 all single particle 
states are metallic, therefore this region is not of interest. 
However for A > 2 single particle states are localized. 
Here we find a large region of metallic bound states for 
3 < U < 15 and A < 3. Towards larger values of A the 
existence region breaks up into two main tongues, which 
we observed also in Figj2] Note also two separate tongue 
structures at around U ~ 1.5 and U ~ 15 which stretch 
up to A w 2.4. 

The observed correlated metallic bound states have to 
form bands with a continuous spectrum. The width of 
these bands will determine the largest group velocity of 
spreading fronts as in Fig(3j We found that these veloc- 
ities are in general depending on the control parameter 
values. In the following we will discuss a possible mech- 
anism for the observed effect of appearance of correlated 
metallic bound states in the main existence region of the 
correlated metal phase in Figj4]for A = 2.5. Let us con- 



FIG. 4: (Color online) The phase diagram of the TIP 
in a quasiperiodic potential as suggested by the long term 
wave packet evolution (the color code maps its second mo- 
ment square root, see text for further details). In light areas 
metallic eigenstates exist, in dark - localized eigenstates only. 
Dashed horizontal line marks the MIT for the noninteracting 
case. Below it the system is a metal. For nonzero interactions 
a new phase of a correlated metal is formed in the midst of 
the insulator region. 



sider the Fock space of two noninteracting particle eigen- 
functions 

M = J—j J2 A t A fKbt\o)- (5) 

where /i < v are the indices of the single particle eigen- 
states, sorted according to their position along the chain. 
For nonzero interactions we expand an eigenstate as 
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and the evolution equation for the coefficients (j)^ v reads 

The renormalized energies E^ v and the overlap integrals 
Ip^/j// responsible for hopping between the modes |/i, v) 
and |//,z/), are obtained from 

£^ = {\» + K) + UI%, 7^ = _^_ E . ( ^)2 

V.^A^V'. (8) 
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We are interested only in bound states where the two 
particles are within a localization length distance from 
each other, since these are observed to yield a transition 
to a correlated metal. The renormalization in each mini- 
band depends on the strength of the overlap integrals I® u . 
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We compute them numerically and find that the average 
overlap integrals 1° ~ 0.5 from TP1 and 1° « 0.3 from 
TP3. Note that the values for 1° for states from TP2 and 
TP4 are much smaller, since the Fock states are made of 
products of different single particle states. TP5 yields 
again large values of 1° but is irrelevant for reasons given 
below. We also obtain that the average overlap integrals 
f^iv^'v' ~ 0.1. Let us consider only states from TP1. The 
increase in U leads to a broadening of TP1 width, since 
some energies get strongly renormalized and some less. 
In the Fock space we therefore observe an increase of an 
effective potential strength (which is similar to A for the 
noninteracting case) as 0.5U. At the same time the differ- 
ent Fock states from TP1 increase their overlap (which 
is similar to the hopping strength of a single particle) 
as 0.1/7. Therefore the increase of the potential strength 
wins, and these states do not cross over into a delocalized 
regime, when no further Fock states are considered. The 
same is essentially true for all TIP minibands. However, 
at some value of U some renormalized states from TP1 
will resonate with weakly renormalized states from TP2. 
This group of states is characterized by a zero potential 
strength (since they are resonant) and any finite overlap 
will therefore lead to a complete derealization. The ex- 
pected value of U follows from the distance between the 
minibands which is around 2, and with 1° = 0.5 we pre- 
dict the derealization to start around U = 4 as observed 
in the numerics. This is the rough location of the left 
large tongue in Figj4j 

Using the same reasoning we predict that a resonant 
mixing of renormalized states from TP3 with ones from 
TP4 is expected at around U = 6.7, as follows from the 
miniband distance around 2 and the TP3 value 1° = 0.3. 
Again this is indeed the observed location of the right 
large tongue in Figj4| It is a challenging task to extend 
the above arguments to the whole phase diagram in Figj4] 

The novel state of a correlated metal formed from two 
interacting particles should be easily measured using in- 
teracting pairs of ultracold Rb atoms in optical lattices, 
which were shown to bind despite their repulsive interac- 
tion [16]. 

To summarize, we observed a nonperturbative der- 
ealization of two interacting particles in a quasiperiodic 
potential deep in the insulator phase of the noninter- 
acting problem. The corresponding correlated metallic 
bound states keep both interacting particles at distance 



less than the single particle localization length. This 
happens because the interacting particles may redis- 
tribute their total energy into different Fock states 
which are coupled due to the interaction. We gave 
estimates for the appearance of the correlated metallic 
phase in the parameter space, which agree well with 
the numerical data. It is a challenging task to extend 
the theory to the whole parameter space, to make 
it quantitative, and to explore the effect of further 
increase of particle number. At finite particle density we 
expect to observe a many body correlated metallic phase. 
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